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Low-complexity Non-coherent Signal Detection for 
Nano-Scale Molecular Communications 


Bin Li\ Mengwei Sun\ Siyi Wang^, Weisi Guo^, Chenglin Zhao^ 


Abstract 

Nano-scale molecular communication is a viable way of exchanging information between nano-machines. 
In this letter, a low-complexity and non-coherent signal detection technique is proposed to mitigate the inter- 
symbol-interference (ISI) and additive noise. In contrast to existing coherent detection methods of high complexity, 
the proposed non-coherent signal detector is more practical when the channel conditions are hard to acquire 
accurately or hidden from the receiver. The proposed scheme employs the concentration difference to detect the 
ISI corrupted signals and we demonstrate that it can suppress the ISI effectively. The concentration difference is a 
stable characteristic, irrespective of the diffusion channel conditions. In terms of complexity, by excluding matrix 
operations or likelihood calculations, the new detection scheme is particularly suitable for nano-scale molecular 
communication systems with a small energy budget or limited computation resource. 

Index Terms 

molecule communications, non-coherent detector, inter-symbol interference, diffusion channel, low-complexity, 
energy efficient 


I. Introduction 


M olecular communications transmit information through the diffusion of chemical molecules [l]-[3]. 

Inspired hy nature, the diffusion channel has advantages in both transmit energy efficiency and propagation 
loss over conventional electromagnetic (EM) wave channels, e.g. the total energy hitting the receiver suffers a lower 
rate of distance dependent loss [4], [5]. Perhaps this is one of the reasons why molecular communication systems are 
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abundant in nature, both at the nano- and at the macro-scale [6]-[9]. For macro-scale applications (e.g. moths and 
crustaceans [7]), the molecular diffusion is often assisted by ambient flows, significantly accelerating the diffusion 
process. Furthermore, the data carried by diffusion channels are delay-tolerant. For nano-scale applications, the 
distance is small and the time delay remains within a few micro-seconds. As for the area of nano-medicine, most 
existing nano-machines lack the wireless communication capability to work effectively towards a common goal, 
e.g., targeted drug delivery or surgery. The aforementioned example of applying efficient micro- or nano-scale 
communications to nano-devices creates a new paradigm to future medical applications [10]. 

Recently, a simple practical testbed has been developed in [11], one that can send generic text information over 
a few meters of unbounded free-space. The implementation consists of a transmitter that emits a short pulse of 
molecules bearing the information data using binary concentration shift keying (B-CSK) modulation, which was 
firstly introduced by M. S. Kuran et. al. [12]. The channel is a random walk diffusion channel with the induced drift. 
The receiver consists of a chemical detector, which wifi sense the concentration of molecules at a specific sample 
time and demodulate the signal. Since the diffusion channel impulse response (CIR) has a long tail, molecules from 
the previous emissions will inevitably intermix with the response of the current pulse, resulting in the inter-symbol 
interference (ISI). In some cases, this can result in aggregated ISI, when the chemical sensor cannot absorb incoming 
molecules [13]. As ISI is the dominant source of errors, one of the major challenges is in designing signal processing 
methods to combat ISI. Recently, the maximum a posteriori (MAP) and minimum mean square error (MMSE) 
schemes have been proposed in [14]. By acquiring the channel response, a decision-feedback equaliser (DFE) is 
designed. Similarly, transmitter side chemical pulse shaping methods proposed in [15] also demand knowledge 
of channel parameters. Such coherent methods could mitigate the ISI, but requiring accurate channel response 
estimation. 

In this paper, we propose a novel non-coherent detection scheme for low-complexity molecular communication 
receivers. In sharp contrast to existing coherent schemes that require the accurate channel state information (CSI), 
the proposed scheme is premised essentially on the difference of the integrated concentration between two adjacent 
bins. The difference in practice is a stable characteristic, irrespective of the diffusion channel conditions. Despite 
its simple implementation, i.e., excluding any matrix operations and likelihood calculations that are indispensable 
in coherent schemes, the non-coherent scheme may also eliminate the ISI effectively. Numerical simulations are 
used to validate our presented scheme. The results demonstrated that the proposed non-coherent scheme may even 
be comparable to the coherent MAP and MMSE algorithms in good operating conditions, but does so with a much 
lower computational burden and requiring no channel information, which is hence attractive to practical molecular 
communication system integration. 



3 


Transmitter 

At) ^ 

Diffusion 



Signal 

shaping 

- W 

Channel 



7 ^ 

Detector 


Fig. 1. A physical-layer schematic structure of molecules communication. 

II. System Model 

In this section, we will first elaborate on the molecular signal modulation and channel model. This is then 
followed by a short introduction on some existing coherent detection schemes, e.g., MAP and MMSE. 


A. Single Model 

A typical physical model of molecular communications is illustrated in Fig. 8. The nano-device transmitter emits 
a short pulse with the duration of Tp, with every T;, interval. The on-off keying (OOK) is used along with a B-CSK 
modulation scheme [16], i.e., the emission is enabled if the binary symbol G A. = {0,1} is 1, /c = 0,1,..., oo. 
Thus, the emitting pulse sequence is: 

= f; Ci X red - trA (1) 

k=Q \ J 

Without loss of generality, in (1) a rectangular pulse shaper rect (t) is adopted, i.e.. 


rect 


t - Tp/2 


'A, -Tp/2 < t < Tp/2, 

,0, otherwise. 


( 2 ) 


where A represents the chemical concentration. After propagated via a diffusion channel h{t) and contaminated by 
other imperfections, the received concentration is: 


“ /1 — T /2 \ 

s{t) = A'^ak- rect I —— ^- kTb | (g) h{t) + w{t) 

k=o V J 

OO 

= A ^ akvit - kTb) + w{t), 

k=0 


( 3 ) 


where the notation (g) denotes the convolution operation; y(t — kTb) = rect[(f — Tp/2)/Tp — kTb] ® h{t) represents 
the equivalent channel response between the binary information source {uk] and the receiver; w{t) is the additive 
noise aroused by the imperfect counting process (when the low complexity chemical sensors with low sensitivity 
are considered), which is assumed to be an independent and identically distributed (i.i.d) zero-mean white Gaussian 
noise [17], with a variance of cr^, i.e., w{t) ~ AA(0, cr^). 














4 


We consider a 3-D environment where the movements of molecules are characterised hy the free diffusion (or 
Brownian motion) [18]. The well-estahlished impulse response of such diffusion channels, hy further considering 
the absorbing effect from a molecule receptor, will be given by [19], [20]: 


h{t) 


mR{d -R) \ {d- Rf 

(i(47rf/)3/2 ^ 


( 4 ) 


where m is the total number of emitted molecules; R is the radius of the absorbing spherical receiver; / is the 
diffusion coefficient of specific medium; d is the Euclidean distance between transmitter and receiver [21]. As 
mentioned previously, unlike classical wireless channels, the diffusion CIR y{t) will exhibit a very long trail given 
the specific d and R. 

For convenience, it is assumed the synchronisation of molecular communication links has been accomplished [22]. 
The nano-receiver will sample the molecules concentration inside the reception space at a Nyquist rate R= l/T^. 
The sampled discrete signal will be: 

OO 

rk= ^rnyk-m + Wk, (5) 

m=0 

where = r{kTb),yk-m = y{kTi, — mTb),Wk = w{kTb). For /c — m < 0, we will have yk-m = 0. If the infinite 
ISI is truncated, the received signal may be expressed as: 

fc-i 

Tk = Okyo + X] Oimyk-l-m + Wk, (6) 

m=k—I 

where the second term accounts for the ISI coming from the slow decay of previous CIRs. In [14], the ISI is 
combined with the additive Gaussian noise, which are then formulated as the non-stationary and signal dependent 
noise. Here, I is defined as the truncating memory length, which specifies the interfering span of ISI from the 
previous transmitted pulses. Note that, whilst the value of I is unbounded in practice, a finite value is used in 
this work for simulation and realistic algorithm purposes. For example, in [14], a finite memory is assumed in the 
sequential MAP and MMSE schemes. 


B. Coherent Detectors 

1) MAP: The received vector is conveniently expressed as r = Ya -|- w, where ri^xi — [t’kj'r'k-i, ■ ■ ■, Vk-K+iY', 
^{K+i)xi — Wk,ak-i, ■ ■ ■ ,o,k-K-i+i]'^', ^Kx{K+i) is the circulant channel matrix constructed from y/xi, and 
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Wifxi is additive noise vector. The channel matrix kx{k+i) is then expressed as: 

2/0 2/1 • • • yi 0 • • • 0 

Y _ 0 2/0 2/1 • • • 2 // 0 0 

0 • • • 0 2/0 2/1 ••• 2/7 


For a coherent MAP scheme, the accurate estimation of unknown CSI, i.e. y/xi, should he acquired firstly. The 
estimation of unknown symbols, relying on the posterior density, is derived from: 

aMAP = argmaxp(a|r, Y), 

K+I K+I ( 7 ) 

= argmax n piak\ai:k-i) p{rk\ri:k-i,Oii:k), 

aeA^+^ fc=i k=i 

where the likelihood density p{rk\ri:k-i,ai:k) follows a Gaussian distribution [14]. To reduce the complexity, a 
sequential estimation can be used in practice [23]. For the i.i.d information source ak G A with the equal priors 
for “0” and “1”, the MAP scheme is equivalent to a maximum likelihood (ML) method. 

2) MMSE: Another widely used sub-optimal detector is based on the MMSE criterion [23], which minimises 
the covariance matrix of estimation errors, i.e.. 


aMMSE = argminE 
aeA^+^ 


(a 


a) (a 



( 8 ) 


Based on the linear Gaussian model in (6), the MMSE estimation of the unknown information is derived via: 


aMMSE = IE(a|r), 


( 9 ) 


= E(a) + (A-i + Y^A-iY)-iY^A-i[r - YE(a)], 

where E(.) represents the statistical expectation; A^ is an K x K diagonal matrix with the elements of its principal 
diagonal are all Note from (9) that, an accurate CIR is also required by the MMSE scheme. 


III. Non-coherent and Blind Signal Detections 

It is recognised that, in sharp contrast to wireless communications, there are two inherent challenges in molecular 
communications. Eirstly, the CSI cannot be easily estimated without a pilot channel, especially in the presence of 
channel disturbances due to air flow and temperature variations. Even in the presence of a pilot channel, the 
coherence period of the channel is usually small, whilst the channel delay is large, making the precise estimation a 
rough work. Secondly, molecular communication is targeted towards nano-scale systems, which naturally demand 
the low-computational complexity architecture to reduce hardware costs and energy expenditure. The effectiveness 
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of the coherent MAP or MMSE algorithms in eliminating ISI largely depends on both accurate CSI and complicated 
matrix or polynomial (likelihood) operations. Based on the above considerations, those coherent detectors will be 
less attractive in the context of low-complexity and low-power nano-machines. 


A. Concentration Difference Estimator 


As far as we are aware of, for the first time a low-complexity non-coherent receiver is proposed for molecular 
communications by this work. The designed scheme essentially involve 3 steps. 

1) Pre-smoothing: In order to suppress the additive noise, firstly a pre-smoothing (or pre-filtering) process will 
be used. To do so, the received molecular concentration is over-sampled around each symbol concentration (e.g., 
kTb) with a higher rate l/Tg, as shown by Fig. 2. Then, a moving average process with a length of 2L -|- 1 is 
applied. Here, we assume LTg = Tb/2. The smoothed signal is given by: 


h,i= 


1 


lo=-L 


2L + 1 


X rk,i+io- 


( 10 ) 


Here, denotes the Ith sub-sample around the A;th symbol {—L < I < L). In practice, L zeros will firstly be 
padded to both the front and the end of received signal sequence, i.e., f = [0,..., 0, ..., ri..., rK,-L, 

rK-L+i, ■ ■ • ,?’iV,L,0, ... ,0]'^. 

2 ) Concentration Difference: In the new non-coherent scheme, the difference between concentration integrations 
is proposed to demodulate information, which is defined as: 

L L 

7fc - XI - X ( 11 ) 

l=-L l=-L 

In the following, we will show that, for a typical transmission rate Tb = Ntm (e.g. N ^ 2) where tm oc d? 
denotes the interval between the maximum concentration from the beginning of CIR, the suggested concentration 
difference (i.e. 7 ?;) may effectively suppress the ISI and achieve the promising detection performance, even if it 
can not mitigate it completely. 
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The molecules residue concentration at the time k, generated from the {k — fco)* previous transmission, will he 
approximated hy: 

Cfc-fco (^) 0 ~ Oik-ko X y\t ~ {k — ko)Th + LTg] + bk-ko x {L 1) + o{{L + 1)“^}, 

( 12 ) 

— Ck-ko + bk-ko X (L + (), — L ^ I < L. 

It is noted that, from the expression of channel response in eq. (4) and the illustration in Fig. 9, the ISI falling 
into subsequent time is not linearly varied with the index k, leading to great difficulties in the algorithm designing 
and analysis. From Eq. (12), it is shown that, fortunately, the ISI from previous molecular concentrations could 
he approximated hy the first-order Taylor series expansion (TSE). Although it may underestimate ISI hy dropping 
the high-order series, for the channel response y{n) with a long trail, the above approximation is sufficient for the 
practice algorithm developments and analysis. With the above Eq. (12), now the nonlinear interference curve will 
be replaced by a linear function of a negative slope, and the ISI from the {k — /co)th time is then approximated 
via the area under this line, i.e., ISI ~ J2i"=-L^k-ko{k,l) ~ ^Si + AS 2 with ak-kg £ {0,1}. Corresponding to 
the first and second terms in eq. (12), the ISI (i.e. a trapezoid area ) can be broken into a square area A5i and the 
other triangle area AS 2 , as in Eig. 9. It is easy to see the constant Ck-ko (1^) accounts for the intercept (or 

the width of a shadow rectangular Si), while bk-kg — Ofc-fco x ^^^1 ,, , in the second term denotes the 


t —(fc— kg')Tj,-\-LTs 


in the second term denotes the 


linear slope. 


Based on Eq. (12), during the /cth time bin, the accumulated concentrations can be further written to: 


X] = X - LTs) + X ek-kg{k,l) 


~ X ~ ~ LTg) + 2L X Ck-kg + 2 X kk-kgL^: 

l=—L ko=l ko=l 

L I 

= otky{t — kTi) — LTg) + 2L Ck-kg I-o{IL‘^). 


useful signal concentration 


ISI concentration 


Here, we omit the additive noise for convenience. The third term may be dropped as the slope is usually very 
small, i.e. bk-kg —^ 0 when the ko is large. The concentration difference 7^+1 is now reformatted by: 

L L 

7k+l = X ^k+l,l — X ^k,l, 

l=-L l=-L ^ 

y X ^k+iy [t- {k + l)Tb - ITs] - X ^ky{t - kTb - ITs), 


where (a) follows the fact that: 


2L X Ck+l-kg Ck-k 


(15) 
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Fig. 3. The illustration of an instantaneous realisation of a molecular diffusion channel y(t), as seen the top sub-figure, and an example of 
received signals r(t) in the below sub-figure. 

From the relationship of eq. (14), it is noted that the suggested metric 7 fc-i-i involves only the differences of the 
useful signals. In other words, the additive ISI has been eliminated basically. A physical explanation to eq. (14) is 
that, as suggested from Fig. 9, the ISI integration fallen into the time bin [kTb — LTg, kTh + LTg] from the (k — koY^ 
molecules emission will be approximated by the combined area of the trapezoid, i.e., + AS' 2 . Similarly, the ISI 

integration fallen to [{k + 1) x T;, — LTg, (fc + 1) x Tf, + LTg] may be approached via S 3 + AS' 4 . It is found that 
the rectangular area Si is almost equal to S 2 and, therefore, we may have Si + AS 2 — (S 3 + AS 4 ) ~ Si — S 2 cr; 0. 

3) Post-smoothing: A post-smoothing process will be used to further reduce additive noise, and a similar moving 
average process may be utilised. Finally, a threshold will be specified to determine whether there contains sufficient 
energy in the current time index k, indicating whether the new molecule concentration appears. In practice, the 
threshold is simply set to 0 . 

B. Complexity Analysis 

From the above elaborations, it is seen that the summation operations of our non-coherent scheme is about 0{L) 
in complexity, i.e., requiring no multiplication. For the traditional MAP scheme, the numbers of multiplication are 
usually measured by 0{2^9), where 9 accounts for the complexity of evaluating the likelihood function [24]. For 
the linear MMSE scheme, the complexity is For the typical configuration of L = 5 and I > 20, we may 

easily note that the proposed non-coherent scheme superior in its low implementation complexity. 
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IV. Numerical Experiments 

In the simulation setup, the maximum concentration of the channel response y{t) is located at tm = as 

seen from the Fig. lO-(a). The duration between two binary symbols is Tf, = 20rs. As shown in Fig. lO-(b), after 
diffused propagation, serious ISI will occur at the receiver. In the analysis, we consider ISI with a relatively long 
interfering length of / = 30. The moving average length is set to L = Each bit error ratio (BER) curve is 
derived numerically based on 10® symbols. Similar to [14], we treat ISI as the receiving self-interference. Given a 
total length of information bits K^, the SNR in the presence of additive counting noise is alternatively defined as: 


SNR = 10 logio 


E 


\ 

,fc=i / 




(16) 


As suggested by [14], there also contains the ISI in the received signal r{k). Another metric, known as the signal 
to noise and interference ratio (SINR), will be also investigated, as given by: 


SINR = lOlogio < 


(Kt, L 

E E yit-m- iTsf 

\k=ll=-L , 


I f: E [rit - m - ITs) - y{t - kn - ITsW 

k=l l=-L 


(17) 


A. Performance Analysis 

From the Fig. 11, we may observe that, given a relatively high data rate, e.g. Tj, = 20Ts = lA3tm, the non¬ 
coherent algorithm will obtain the promising detection performance. It is found that, by suppressing the ISI to the 
minimum, the sequential MAP scheme will achieve the best performance in the high SNR regime. Although the 
proposed non-coherent detector has a very low-complexity implementation, the new scheme is able to achieve a 
superior performance in comparison with the coherent MMSE scheme, but cannot achieve a better performance 
than the coherent MAP scheme. More importantly, as a convenience the CSI can be excluded from our blindly 
non-coherent detection scheme, unlike classical coherent MAP or MMSE schemes. 

In practice, the symbol duration Tf, may have significant influences on the detection performance. Firstly, we note 
from Fig. 9 that, in the low SNRs region (e.g., <16 dB with = 20Ts), the non-coherent scheme will surpass the 
two coherent schemes. In high SNRs region, however the coherent MAP scheme will be superior to the non-coherent 
detector. Primarily this is because the designed non-coherent scheme is relatively robust to additive noises (e.g. 
with a moving average process), while the coherent schemes are more sensitive to such contaminations. Secondly, 
it is observed that, as the symbol duration decreases (or the transmission rate increases), the performance of all 

'Note that, here we only focus on the processing of discrete samples. While the sampling time T, may depend on different applications, 
the signal processing scheme is independent of any specific Ts- 
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Fig. 4. BER performance of coherent detection schemes (MAP and MMSE) and the non-coherent scheme symbol interval. Note that, we 

r '^h 

assume L = . 


three methods will he degraded. Taking L = 7 {Tb = lATg = tm) as an example, the detection performance of the 
non-coherent scheme will further he worse than both MAP and MMSE techniques. This is mainly attrihuted to the 
fact that, with a small L, the ISI of two successive time-hins (especially when /cq = 1) could not he approximated 
hy its first-order TSE anymore. In such a case, the concentration difference may not mitigate the ISI as effectively 
as set out in Eq. (14), leading to the serious performance deterioration. Eor coherent MAP and MMSE schemes, 
the residual errors will also he more apparent in the presence of more severe ISI. 

Eor some specific molecular applications calling for the low implementation complexity and low operational power 
expenditure, the memory length of ISI in algorithm processing should he restrained when implementing MAP or 
MMSE algorithms. If the processing ISI length is restricted to I = 10, from Eig. 12 the proposed non-coherent 
scheme will surpass both MAP and MMSE detectors, resulting in a superior ISI mitigation performance. 

Eurthermore, as the ISI is treated as the receiving interference, we may investigate the BER performance of 
different detection schemes in the context of SINRs. As seen from Eig. 10, the channel response y{n) is assumed to 
exhibit a long trail, the self-interference from previous transmissions, therefore, will dominate the molecule receiver. 
In addition, we have considered another Gaussian noise to characterise the sensing (or counting) errors of a chemical 
receptor. It is found from Eig. 13 that, given a long-trail ISI response (e.g. I = 30) and the equal probability of Is 
and Os in an OOK line-code, the SINR will converge to a constant value (e.g. —3 dB for the considered symbol 
period) as the additive error tends to be 0, indicating the ISI will always surpass useful signals, which may remain 
quite different in comparison to wireless communications. Erom the simulation result, we observed that, in realistic 
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Fig. 5. BER performance of coherent detection schemes (MAP and MMSE) and the noncoherent scheme under various I. 
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Fig. 6. BER performance of various detection schemes under different SINR. 


molecular communications, the proposed non-coherent detection scheme will become comparable to both MAP and 
MMSE as the SINR approaching —3 dB (or the variance of additive noises approaching 0). 


B. Comparative Analysis 

Note that, the modulation scheme investigated in the above analysis can be viewed as a special case of the binary- 
concentration shift keying (B-CSK) [12], i.e., no molecule will be transmitted in the case of “0”. As suggested by 
ref. [12], a one-step sequential interference cancellation (SIC) scheme would be adopted to mitigate ISI from the 
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Fig. 7. BER performance of different modulation schemes. 

last time-slot. On this basis, the information capacity has also been derived based on a simplified channel model, 
provided the interference to the time k from previous slots k — ko (ko >2) could be ignored. These previous works 
in [12], [25] may be of great importance to the low data-rate molecular communications (i.e., the symbol duration 
is greatly larger than tm)- As far as some scenarios of relatively high data-rate are concerned, however, the ISI 
from multiple previous slots will occur at time k, as seen from Fig. 10, which would make the one-step SIC 
detection infeasible. Given the ISI channel of / = 30, the numerically derived BER curve are plotted in Fig. 14. 
It is seen that, for one thing, the one-step SIC scheme fails to deal with a long ISI channel, and the BER value 
can only decrease to a floor of 7.4 xl0“^. For another, even if the SIC method is generalised to multiple steps, 
the BER performance will be converged to a floor of 6.3 xl0“'^ if the SIC span surpasses 5. It is considered that 
the BER floors are mainly aroused by the error propagation of SIC. That is, once a detection error occurs, it will 
affect the detection of subsequent time slots. Besides, the threshold configuration after multiple-steps SIC requires 
full information of CIR, making it less attractive in low-complexity applications. From the numerical results, the 
proposed non-coherent detection scheme will be inferior to the multiple-steps SIC in low SNRs. By mitigating the 
BER floor effectively, however the new non-coherent scheme will be more competitive in high SNRs. 

For the comprehensive analysis, we further derive numerically the BER performance of another modulation 
scheme, i.e., binary molecule shift keying (B-MoSK) [12], [25], in which two distinct molecules are utilised to 
convey two different information bits (or messages). It is seen that the interference from the previous time slot (i.e. 
k — 1) would be fully mitigated under a simplified channel model (i.e. with the Ist-order ISI). For the ideal case, the 
B-MoSK modulation scheme will achieve the promising detection performance by mitigating ISI via modulation 
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techniques, as demonstrated by Fig. 14. When the long ISl channel is considered (i.e. corresponding to a small 
symbol duration), however the interference from previous time (A: — 1 — ko) (ko > 1) may appear at the current 
time k. Even if a similar SIC technique can be applied, the BER performance of the B-MoSK modulation will be 
degraded then. Moreover, with the integrated two different chemical receptors, the implementation of a B-MoSK 
receiver may easily be impractical to certain low-complexity scenarios. Meanwhile, note that if these two molecule 
receptors are not completely separable, then the adjacent information corruption will further occur. In contrast, 
the non-coherent scheme, which allows a compromise between the complexity and performance, will be of great 
promise to the low-complexity and robust molecular communications. 

V. Conclusions 

A low-complexity and non-coherent signal detector is proposed to mitigate ISl in the diffusion-based nano-scale 
molecular communication systems. This new scheme is premised on the difference of molecules concentrations, 
which is a stable characteristic irrespective of channel conditions. Therefore, this method excludes the need for 
accurate and timely channel estimations. Both the theoretical analysis and experimental simulations demonstrate 
that, provided a moderate transmission rate (e.g., = lAtm) and a long ISl channel, the proposed non-coherent 

scheme is able to achieve the comparable error rate performance compared with coherent MMSE or MAP schemes. 
The proposed non-coherent detector is of great promise to low-complexity and low-rate molecular communications, 
as it completely excludes the matrix operations and likelihood calculations and thereby has the low computational 
and implementation complexity. Euture work includes its extensions to extremely high rates (Tf, < tm) and the 
hardware implementation. 
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